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Abstract. In this paper, I explore various transport processes that have 
a large impact of the distribution of elements inside stars and thus, on 
stellar evolution. A heuristic description of the physics behind equa- 
tions is provided, and key references are given. Finally, for each pro- 
cess, I will briefly review (some) important results as well as discuss 
directions for future work. 



1 Atomic diffusion 



> 
On 

'1.1 Introduction 



OO ' Atomic diffusion refers to atomic scale physical processes that cause the strati- 
fication of various species and the smoothing of chemical discontinuities. These 
processes occur when multicomponent plasmas are subject to physical conditions 
that vary within the medium (such as pressure, temperature or density) or are 
subject to external forces that vary from species to species and are related to the 
Brownian motions of particles. Stellar interiors are a good laboratory for the study 
!— J ' of these processes. 

As atomic diffusion is not the main subject of this paper, this section will briefly 
outline general ideas to keep in mind when diffusion^] is mentioned. Calculations 
are presented for a small concentration of particles i diffusing in protons p. This 
is called the trace element approximation. For a more complete description of 
diffusion in a general context see Burgers (1969) or Chapman & Cowling (1970). 
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1 het us note that atomic diffusion is generally referred to simply as diffusion. 
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1.2 Diffusion and external forces 

We first consider a gas in which physical parameters are constant and some external 
force .Fextj is exerted on particles i. The total force on i is then 

r i = -k B T?^+r estl (i.i) 

or 

where Pi represents the partial pressure on i, and is related to the total pressure 
P and the concentration Cj by 

Ci = P ' with Pi « P p . (1.2) 

*p T" *i 

Equilibrium is obtained when Ti = 0. Since diffusion is a slow process that should 
bring composition gradients to their equilibrium values, it is reasonable to assume 
that the diffusion velocity Vi p of i with respect to protons is proportional to the 
difference between actual and equilibrium values (Eddington 1930) and thus 



d In Ci T c: 



dr k B T 

and the change of c$ during this evolution is given b>H 



(1.3) 



The atomic diffusion coefficient is D at = ^Avt with A the particle mean free path 
and Vt the thermal velocity (see e.g. Aller & Chapman 1960). In the absence 
of an external force, this equation shows the tendency of microscopic motions to 
eliminate abundance gradients. 



1.3 Gravitational settling and radiative forces 

Let us now look at two specific external forces encountered in stellar interiors. 
First, we have gravity. When the mean molecular mass of the stellar plasma /i is 
different from the mass of i (rrij), the differential force on i is 

•^exti = M ~ m i9 = (M - m i)9 (1-5) 

where we used again the perfect gas law, and assumed hydrostatic equilibrium, g is 
the local gravitational acceleration. In stars, this equation implies that elements 
heavier than hydrogen must "fall" towards the star's center, and this is called 
gravitational settling. 

The other force that has a large impact on element stratification is the transfer 
of momentum caused by photon absorption. In the region where radiative transfer 



2 Here, modifications of the structure are neglected. 



S. Talon: Transport processes in stars 



3 




Fig. 1. The absorption of photons of different energies by various atomic species causes 
an acceleration that opposes gravity. When absorption is more probable, and the number 
of absorbers is small, this may produce their radiative levitation. 



can be described using the diffusion approximation, the star's thermal structure 
is calculated using the Rosseland mean opacity kr defined as 

i r dB{ ,,T) d ^ r- L o_Bpr) dv (16) 



kr Jo dT J k v dT 

where B(u,T) is the Planck function and k„ the monochromatic opacity. The 
radiative luminosity L lad is then related to the temperature gradient by 

L rad = 647T(7 fl r 2 T 3 dT 

3 K R p dr 

where or is the Stcfan-Boltzmann constant. For a given element i, the absorbed 
momentum is proportional to the fraction /j of the total flux that is absorbecH by 
that type of particle. We can now compute the radiative acceleration of particle i 

as 

1 i rad h 



where 



9 ^-a^— Kr Y, (L8) 



K ' J, -V(u)du with / V{u)du = 1 (1.9) 



(I 



where u = hv/ksT and V{u) = j^j^r^ijz is the normalized radiation flux associ- 
ated with the temperature derivative of the Planck function (this was derived first 
by Michaud et al. 1976). The calculation of radiative forces thus requires detailed 
spectra for all species and all ionization states. 



3 This is in fact an approximation. The complete treatment involves momentum redistribution 
between ions and electrons (see Richer et al. 1998 for details). 
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Fig. 2. When temperature is not homogeneous, the particles mean free path varies. In 
this cartoon, the probability of crossing the line to the left is different from the proba- 
bility of crossing it to the right if the abundance is homogeneous. The achievement of 
equilibrium requires that there be a lower density of big particles on the left than on the 
right. 



For a given trace element, when g ra( j > g, the element will be "supported" by 
radiation, and will move upwards toward regions of force neutrality (g ra d — <?); this 
neutrality condition may also result from absorption saturation^. In the opposite 
case (grad < g) , elements sink. 

1.4 Diffusion and temperature gradients 

Temperature gradients are another source of element stratification. In a binary 
fluid (with all other physical parameters constant), the particles mean free path 
depend on temperature through their collision cross section (roughly as a oc T~ 2 
for ions, while it does not depend on T for neutral atoms) and through their 
thermal velocity (vt oc T 1 / 2 ). 

The efficiency of this process can be estimated by considering the net flux of 
ions through an equipotential (Vauclair, 1983) 



1 , , el, . 1 a 



fi ir+A "«lr+A | fi u + lr- A "»lr-A 

where rij is the number density and assuming that on average 1 /6 of ions move in 
the direction of the equipotential. If is fixed, this is simply 

v, . * »£f> a -inJJf (,.„, 

3cr or 2 or 



4 As atoms migrate upwards and Xi rises, they have to share the momentum that photons 
give away, potentially leading to a saturation of <? rad which limits further growth of Xi. 
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2.50T5.30D10K-3 255 My, T = 9784 K . . 

2.50T5.30D10K-4 253 My, T* = 9752 K — — 

2.50T5.30D1M-4 253 My, T* = 9830 K .- — 

2.50T5.20D1M-4 251 My, T = 9808 K 

2.30R1K-2 250 My, T = 9458 K . 
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Fig. 3. Surface abundances in Sirius under the effect of atomic diffusion including ra- 
diative forces and some ad-hoc turbulence (various lines) compared with observations by 
various authors (symbols). (Figure 18 of Richer et al. 2000.) 



for ions. A realistic calculation would give a leading factor of 2.65(Zi/Z p ) 2 instead 
of the 2 given here. In more general conditions, even the sign in Eq. may be 

different as it depends on the masses of the colliding particles. This is discussed 
at length for example in Michaud (1993). 

1.5 Some applications of diffusion theory 

Diffusion theory has been applied to a wide variety of stars. Its effects are expected 
to be largest in stars which have only shallow surface convection zones. This is the 
case in particular of AmFm stars, which exhibit abundance anomalies that can be 
as large as one order of magnitude for certain elements and thus, it is believed that 
these are only superficial. The Montreal group calculated stellar models including 
all effects of atomic diffusion, including radiative accelerations (Richer, Michaud 
& Turcotte, 2000). One of their results is shown on Fig. [3] When radiative forces 
alone are considered, the predicted anomalies are too large by a factor of ~ 3. In 
the results shown here, a small amount of ad-hoc turbulence has been added to 
obtain a good agreement with the observations. 

Diffusion theory has also been applied to solar models. There, it predicts the 
settling of helium below the convection zone: such settling actually improves the 
agreement with helioseismology as was originally shown by Christensen-Dalsgaard, 
Proffitt & Thompson (1993). The impact of thermal diffusion on solar models has 
been studied with great care by Turcotte et al. (1998) who showed that, for some 
elements, it increases the settling velocity by as much as 50%. 
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Fig. 4. As a fluid element is displaced upwards, it will become denser than its surround- 
ings and will feel a downward restoring force if the Brunt- Vaisala, frequency is positive, 
that is if the temperature gradient is smaller than the local adiabatic gradient. 

Finally, diffusion can also have significant effects on the stellar structure when 
slow changes can accumulate over the stellar lifetime. In old population II stars, 
it is found to reduce the age determination by isochrones in globular clusters by 
about 10% (VandenBerg et al. 2002). 



2 Rotational mixing 

2.1 Introduction 

Descriptions of rotational mixing have evolved a lot in recent years. The need for 
extra mixing in stars is widely recognized, and different rotational histories could 
well explain the variety of stellar behaviors. 

Early descriptions have relied on very simple power laws relating the amount 
of mixing to the star's rotation rate (see e.g. Zahn 1983). In the case of lithium 
destruction in low mass stars for example, this would lead to larger dispersion in 
Li than observed (Schatzman & Baglin 1991). 

Very early on, it was recognized that such scaling laws fail to produce an 
appropriate depth dependent diffusion profile and that consistent models must be 
built that take into account the angular momentum evolution (Endal & Sofia 1978). 
To do so, one must consider all possible transport processes for angular momentum 
and then analyze their effect on the transport of chemicals. 



2.2 Rotation driven instabilities 

In this first section, we will be looking at hydrodynamical instabilities. They are 
presented by considering what happens if one displaces a fluid parcel from its 
original position. 
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Fig. 5. To take into account the stabilizing effect of rotation, we expand a torus by an 
amount £. 



2.2.1 Brunt- Vaisala frequency 

Let us displace upwards a parcel of fluid from its equilibrium position in a radiative 
zone (Fig. 2]). If the displacement over a distance Sz is adiabatic and if pressure 
equilibrium is maintained, the density inside the parcel p pa rcci will be larger than 
the local density piocai by an amount 



Pparccl Plocal . 
— = Sz 

PO 



( dlnT 




dlnT\ 


V dz 


H 

ad 


dz ) 



_ din it Sz ._ _ _ , 

St—T 1 = T7- V ad - V + V M 
dz Hp 

(2.1) 

where Hp — — dz/dlnP is the pressure scale height. We also used the classic 
notation of stellar modeling V = dlnT/dlnP. The parcel then experiences a 
restoring force g (p pa rcoi — PiocaOi giving rise to an oscillation characterized by the 
Brunt- Vaisala frequency N defined as 



N 2 = AT| + Nl = JL (Vad - V + V M ) . (2.2) 
Hp 



The region is locally stable against convection if N is positive 



2.2.2 Solberg-H0iland instability 

The next instability is related to the Coriolis force (Rayleigh 1916, Taylor 1923). 
We consider a fluid in cylindrical rotation. In the absence of gravity, pressure 
gradients balance centrifugal acceleration 

1 d P 

---^- + sft 2 = 0. (2.3) 
p ds 



Let us displace a fluid torus over some distance £ (Fig.JSJ). The net force acting on 
the fluid is proportional to the difference between the torus' angular momentum 
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Fig. 6. In a baroclinic state, equipotentials and isentropic surfaces do not coincide. In a 
displacement from A to B, the local entropy rises, and this state is stable. However, in a 
displacement from A to C, the local entropy decreases, and this characterizes baroclinic 
instabilities. 

(s+£)^torus aiul the equilibrium value (s+£)f^ + £. Since the Coriolis force ensures 
(specific) angular momentum (j = s 2 fl) conservation, it may be written as 

(s + O s as 

valid to first order in £. We use this condition to define the Rayleigh frequency 

Af 2 = l-f( s ^) 2 (2.5) 
as 

which corresponds to a restoring force if positive, that is, if the specific angular 
momentum increases outward. 

If gravity is taken into account, there are now two restoring forces. For equa- 
torial radial displacements, the stability condition then reads 

N 2 + N^> 0. (2.6) 

For general axisymmetric perturbations, the stability criterion is twofold: 

• Condition (|2.6|) must be satisfied; 

• Specific angular momentum must increase from pole to equator. 

The second condition corresponds to applying the Rayleigh criterion (Nq > 0) to 
an isentropic surface; in the limit of no rotation, the first condition corresponds to 
the Schwarzschild criterion. These conditions have been derived by Solberg (1936) 
and H0iland (1941). 

2.2.3 Baroclinicity 

If the star's rotation state is not cylindrical, it can be shown that surfaces of 
constant pressure, density, entropy and gravity do not coincid<H. This is called a 



5 Baroclinicity could also be caused by uneven heating (in a binary system ) and this is the 
case in the Earth's atmosphere. 
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baroclinic state (in opposition to a barotropic state). Let us now imagine displacing 
a fluid element verticalljQ in such a star, going from A to B in Fig. [5] Assuming 
adiabaticity, the fluid will then be denser than its environment and gravity will 
restore it to its original position. However, if the same parcel is moved to C, 
it will now have a lower density, and gravity will make it rise, leading to an 
instability, which is called a baroclinic instability. However, angular momentum 
conservation makes axisymmetric movements towards C impossible (Knobloch & 
Spruit 1982). In the case of non-axisymmetric disturbances, fluid element exchange 
angular momentum (Cowling 1951), allowing instability, which gains energy by 
lowering the center of mass of the corresponding fluid. 

These non-axisymmetric perturbations have been studied in the geophysical 
context using the /3-plane approximatiorQ- Knobloch & Spruit (1982) present a 
summary of these results, and conclude that instability is possible only for fast 
rotation that is if 



2.2.4 GSF instability 

In stellar interiors, thermal diffusivity Kt is much larger than viscosity v which 
controls angular momentum diffusion. This weakens the stabilizing effect of the 
density stratification (Goldreich & Schubert 1967, Fricke 1968). In the case of a 
vanishingly small viscosity, Goldreich & Schubert (1967) concluded that 

... a necessary condition for stability is that the angular momentum 
per unit mass be an increasing function of distance from the rotation 
axis 

which can also be expressed as 



which is similar to the Rayleigh criterion for incompressible, inviscid fluids (see 
e.g. Chandrasekhar 1961). This instability which depends on heat diffusion has 
a lower growth rate, and this changes the nature of the instability (it becomes 
"diffusive"). Indeed, for thermal diffusivity to be efficient, the length scale of the 
perturbation has to be small, and this reduces the efficiency of the instability. 

Fricke (1968) and later Acheson (1978) considered the effect of a finite viscosity 
and presented and a modified version of the Solberg-H0iland criterion for stability 



6 The vertical is defined by the direction of gravity. 

7 This approximation consist in a local study of the equations in a tangent plane touching the 
surface of the sphere at a given latitude. The Coriolis force is then expanded linearly around its 
local value fi = 2f2 cos 6. 




(2.7) 



N&>0 and — = 0, 



(2.8) 




(2.9) 
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together with 

an 2 
dz 



< -£-N*, (2.10) 



permitting a small amount of differential rotation in z before instabilities set in. 
2.2.5 ABCD instability 

If the star is not chemically homogeneous, mean molecular weight gradients also 
contribute to the density stratification. Just like the restoring force associated 
with the temperature stratification is reduced by the thermal diffusivity, this con- 
tribution is also weakened, but this time by molecular diffusivity K^. The Solberg- 
tfoiland criterion is further modified and becomes 

-^N* + ^-Nl + Nl>Q. (2.11) 

Since is rather small (of the order of v), this instability likely does not set in. 
There is however another possibility. When displacing a parcel towards C, it will 
come back to its equilibrium position A due to momentum conservation. However, 
if thermal diffusion is large enough, it will then be cooler than its surroundings, 
and the parcel will sink further in. This produces growing oscillations, a situation 
called over stability. The condition for stability is then (Knobloch & Spruit 1983) 

^- (iV| + Nl) + Nl > 0. (2.12) 

All criteria described so far are linear, and little work has been done on the 
non- linear behavior of the corresponding instabilities, implying that we do not 
know how efficient they are for mixing. 

2.3 Shear instability 

Rayleigh (1880) derived a linear stability criterion against shear, which reads sim- 
ply 

£-» ( , 13) 

everywhere, or in words, there must be no inflection point in the velocity profile 
for it to be stable. In experimental cases however, it is known that turbulence 
can develop even in the absence of such an inflection point (that was shown for 
example by Wendt (1933) in the study of the flow between two cylinders). This 
instability is linked to the presence of finite amplitude fluctuations that can appear 
in the fluid. 

To look into the non-linear regime of interest, we can study the energetics of 
the problem. It can be shown that the minimum energy state of a rotating fluid 
is solid body rotation. Thus, if the star is rotating differentially, by homogenizing 
the velocities, it is possible to extract energy. However, to do so, one must be 
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Fig. 7. To derive the instability condition against shear, let us take a blob at level z with 
a velocity u and move it upward to the level z + Sz. Simultaneously, we take a blob at 
level z + 5z with a velocity u + Su and move it down. When doing so, we also homogenize 
the velocities that become u + 8u/2 for both blobs. We gain 1/ApSu 2 in kinetic energy, 
but lose gApSz in work against gravity. This leads to the Richardson criterion. 

able to overcome the density stratification. By comparing both and assuming 
adiabaticity, one obtains the Richardson instability criterion (Fig. [7]) 



(see e.g. Chandrasekhar 1961). This becomes our non-linear condition for insta- 
bility. 

2.3.1 Dynamical shear instability 

Along isobars, there is no restoring force, and thus, shear is unstable as soon as 
horizontal differential rotation is present. Horizontal shear is thus a dynamical 
instability, acting on a dynamical hence quite short time-scale and leading to a 
large horizontal turbulent viscosity. A few attempts have been made in order 
to evaluate this value (e.g. Richard & Zahn 1999, Maeder 2003, Mathis, Palacios 
&; Zahn 2004), but as of now, there is no reason to prefer one formulation over 
the other. It has strong implications on the transport of elements and we will 
discuss these in § 12.4.41 and 12.61 (this is also discussed at length in Talon, Richard 
& Michaud 2005). 

2.3.2 Secular shear instability 

In the direction of entropy stratification, both the thermal and the mean molec- 
ular weight stratifications hinder the growth of the instability. However, thermal 
diffusion and horizontal shear act as to reduce the stabilizing effect of the stratifi- 
cation, just as in the case of the GSF instability (§ !2.2.4p . We obtain the instability 



N 2 



(2.14) 



Ri = 



(du/dz) 



2 < ffic "t - 7 
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criterion 

(f^) w Hf^>'2<Me) 2 <2 ' 15) 

where T = v£/K T and T p = v£/K M (Talon & Zahn 1997, see also Maeder 1995). 
We normally use i?i cr it = 1/4, but this is open to discussion (see e.g. Canuto 1998) 
This criterion implies that, the smaller the eddy, the more efficient thermal dif- 
fusivity. In fact, there always exists an eddy that is small enough so that the 
instability criterion (|2.15p will be satisfied. Turbulent diffusivity corresponds to 
the largest eddy satisfying (|2 . 1 5[) ; this is a non- linear description, based on energy 
considerations. It is the formulation we shall adopt to describe the vertical shear 
instability. 



2.4 Meridional circulation 

The Eddington- Sweet meridional circulation is related to the thermal imbalance 
that is in general present in rotating stars. This imbalance was first pointed out 
by von Zeipel (1924) when he tried to describe thermal equilibrium in a rotating 
star in solid body rotation. Since such a star is barotropic, equi-potentials and 
iso-therms coincide, and are ellipsoids. When writing the divergence of the heat 
flux 

-V-F = V-(xVT) = -pe, (2.16) 

von Zeipel realized that, while the (horizontal) average of V • F may equilibrate 
the energy generation pe, it is not the case along an equipotential, except if this 
energy generation has a very special behavior, which is not realistic (this divergence 
is shown on the left of Fig. [H in the case e = 0), and this became known as 
von Zeipel's paradox. Vogt (1925) and Eddington (1925) found a simple solution 
to the problem, which consists in considering the advection of entropy S in the 
heat equation 

pTu ■ VS* = V • f%VT^ + pe, (2.17) 

where u is the meridional circulation velocity, which is calculated with the above 
equality. This was done initially for the case of solid body rotation by Sweet (1950) 
and yields a large scale circulation that rises at the pole and sinks at the equator. 

Sweet's solution is however not fully consistent since this large scale circulation 
advects angular momentum as well as entropy. This will modify the rotation 
profile, and hence feedback must be treated. Since this is actually (at least) a 2-D 
effect, not much progress has been made in incorporating it into stellar evolution 
codes until quite recently. 

In a first step, models where built by treating the advection of momentum as 
a purely diffusive process (Endal & Sofia 1976, 1978, Pinsonneault et al. 1989, 
Langer 1991, these codes are still in use). 

The key ingredient to allow the modeling of meridional circulation as a 1- 
D process is the assumption of highly anisotropic turbulence, with the vertical 
turbulent viscosity v v being much smaller than the horizontal turbulent viscosity 
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Fig. 8. The traditional approach to meridional circulation consists of calculating the 
thermal imbalance for a given rotation state (left, here in the case of solid body rotation 
with an oblateness of 10%) and evaluating the advective flux required to counter-balance 
this imbalance (right). 



Vh- This property is related to the vertical stratification (Tassoul & Tassoul 1983, 
Zahn 1992) and permits to assume a state of shellular rotatiord, i.e. il = Q(P). 
Conditions for this approximation to be valid are discussed by Tassoul & Tassoul 
(1983), Zahn (1992) and Charbonneau (1992 - this discussion is based on numerical 
simulations) . 



2.4.1 Heat flux 

Under that assumption, and to first order, it is possible to derive an equation for 
the thermal imbalance in a star with a given (shellular) rotation profile. Then, 
as in the case of solid body rotation, meridional circulation is obtained by stating 
that advection of entropy S must neutralize this thermal imbalance 



pTu-VS = V- xVT +pe-V-F, 



(2.18) 



(this form of the equation also contains a turbulent horizontal heat flux Fh, Maeder 
& Zahn 1998). The vertical component of the asymptotic circulation velocity is 
then 



P 



1 



pgCpT (V ad - V + V„) 



M 



Ea + E„ 



TCpdQ 



S dt 



P? 



(2.19) 



8 Note that, in the case of shellular rotation equi-potentials do not exist, and will be replaced 
by isobars. 
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Fig. 9. Meridional circulation has a different effect on angular momentum and on chem- 
icals. Even when horizontal mixing is strong, the large scale circulation produces a net 
flux of s 2 fi, while there is no net flux of concentration. In fact, the evolution of c is 
possible only in the presence of horizontal inhomogeneities. (From Talon 1997) . 



(see Maeder & Zahn 1998 for the complete expression). The circulation velocity is 
related to a term Eq which depends on the rotation rate, a term Eq which depends 
on the radial differential rotation (9 oc dil 2 /dr) and a term which depends 
on the horizontal variations of the mean molecular weight /i. In a homogeneous 
solid body rotating star, only the first term survives, and one gets Sweet's (1950) 
classical solution characterized by a large circulation cell rising along the pole, 
with a reversed cell at the surface (Gratton 1945, Opik 1951). 

Let us note that if convective overshooting is not present, the solution diverges 
at the convective boundary where V ac i — V = and viscous boundary layers are 
then required (Tassoul & Tassoul 1982). 



2.4.2 Angular momentum transport 



If the only transport processes present are meridional circulation and hydrody- 
namical instabilities, the distribution of angular momentum within the star will 
evolve according to an advection-diffusion equation 



4 w 



J_d_ 

5r 2 dr 



r 4 dr 



pv v r 



dr 



(2.20) 



(see Zahn 1992 for details). Note that the first term of this equation still corre- 
sponds to an advective process since, even if f2 is homogeneous on isobars, s 2 £l is 
not (see also Fig. [9J. This will not be the case for chemicals (see § I2.4.4j) . Here 
we use a vertical turbulent viscosity v v which is attributed to relevant rotational 
instabilities detailed in ^ 12.21 and 12.3.21 Neglecting the star's evolution, this equa- 
tion admits a stationary solution in which advective transport is counterbalanced 
by turbulent transport. If no turbulence is present and if the star is homogeneous, 
circulation will simply stop, resulting in a rotation profile with a core rotating 
about 20% faster than the surface (Urpin et al. 1996, Talon et al. 1997, see also 
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r R 

Fig. 10. Equilibrium rotation profile as given by Eq. (|2.20[) for various surface velocities 
in a 9 Mq main sequence star. (Figure 1 from Talon et al. 1997.) 

Fig. I10p . The equilibrium profile in Q thus depends on the nature and the vigor 
of rotational instabilities present in the star. However, when using realistic values 
for u v , it varies only slightly, and so, meridional circulation is in general much 
weaker than predicted by the simple Eddington-Sweet solution. For any rotation 
profile, the direction of the circulation can be deduced from the requirement of 
reconstructing this equilibrium solution. 

If the star is not chemically homogeneous, meridional circulation will pro- 
duce horizontal variations of mean molecular weight fj, out of the vertical chem- 
ical stratification. These contribute to the heat flux, via the equation of state. 
Equation (|2. 19[) shows that, if horizontal variations of mean molecular weight are 
present, the equilibrium profile described in the previous paragraph exists; how- 
ever, it now describes the function — pt//i (p refers to the horizontal fluctuation 
of fj,). If rotation evolves solely under the action of meridional circulation and tur- 
bulence, this equilibrium requires an increase in differential rotation (compared to 
the homogeneous case) in order to maintain the asymptotic solution. This effect 
has been observed in fully consistent unidimensional calculations including differ- 
ential rotation and mean molecular weight gradients (Talon et al. 1997, Palacios 
et al. 2003). 

If the star is constrained to rotate as a solid bodj@ and does not evolve, the 
build up of horizontal mean molecular weight fluctuations can thus also lead to a 
circulation free state, as originally suggested by Mestel (1953). 

In an evolving stellar model, the required equilibrium profile also evolves (e.g. 
Talon et al. 1997). This leads to the appearance of a "re-adjustment" circulation, 
which is added to the asymptotic circulation. This circulation cannot be avoided, 
even in chemically inhomogeneous stars rotating as solid bodies. 



9 In solid body rotation, no hydrodynamical instabilities can develop. 
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Fig. 11. The combination of advection and horizontal transport leads to an effective 
vertical diffusion. In this cartoon, we are slowly building a discrete version of the error 
function, characteristic of a diffusion process. 



2.4.3 Wind-driven circulation 

If the star's surface is braked via a magnetic torque, as is the case of Pop I stars 
cooler than about 7000 K, the internal distribution of angular momentum is rapidly 
moved away form its equilibrium profile, with large shears developing in the outer 
regions. This is the regime of "wind-driven" circulation, and the model presented 
up to now predicts strong mixing related to such a state. 



2.4.4 Chemical transport 

While momentum transport remains a truly advective process when considering 
the combined effect of meridional circulation and horizontal diffusion, it is possible 
to understand hcuristically that this combination leads to (vertical) diffusion in 
the case of chemicals (Fig.[TTj see also Chaboyer & Zahn 1992; this has also been 
shown in numerical simulations by Charbonneau 1992). This is the case because 
horizontal turbulence continually reduces horizontal chemical inhomogeneities, im- 
peding efficient advection. This effect is felt as soon as 



D h > \rU\ and D h /D v > 1000 



(2.21) 
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(Charbonneau 1992). In the following, we assume that diffusion coefficients are 
equal to the turbulent viscosities i.e. D v = v v and — vh- The advection + 
horizontal diffusion + vertical diffusion equation can then be replaced by a purely 
diffusive equation, with vertical turbulent transport and an effective diffusivity 
depending on horizontal turbulent transport and advective velocity u 



dc 

p at =pc 



1 9 r 2 rr 1 19 

-^7T V P V w c \ + -2 7T 
r z or 1 J r z or 



dc 

r 2 p {D cS + D v ) — 



(2.22) 



The effective diffusivity is D cS = r 2 u 2 / (30A0 (Chaboyer & Zahn 1992); c is the 
nuclear production/destruction rate of the species, and Vi p is the atomic diffusion 
velocity (see §[1]). 



2.5 Some applications of rotational mixing 

This theory of rotational mixing based on an advective formulation of meridional 
circulation and shear instabilities has been applied to a wide variety of stars. In 
the case of massive stars, extensive studies have been conducted by the Geneva 
group, with great success in explaining problems such as: 

• He and N overabundances in O- and early B-type stars and in super-giants; 

• the number ratio of red to blue super-giants; 

• the Wolf-Rayet to O-type stars ratio. 

Recent results regarding massive stars are given in Meynet (2007, this volume). 

In the case of low mass stars, rotational mixing has been applied to the case 
of the lithium dipf 1 "!. In this model, the hot side of the lithium dip corresponds 
to stars that suffer magnetic braking (see § 12.4.31 and 13. 4j) ; this induces large 
internal differential rotation, increasing rotational mixing and thus, surface lithium 
depletion. Rotational mixing reproduces well the hot side of the dip (Talon & 
Charbonnel 1998, Palacios et al. 2003, Fig. [P2f but fails for temperatures below 
T off < 6700 K. 

Rotational mixing has also been applied to the study of the solar rotation 
profile. According to helioseismology, the solar radiative zone is rotating more or 
less as a solid body (see e.g. Chaplin et al. 1999, Couvidat et al. 2003). In that 
case, complete self-consistent calculations predict that radial differential rotation 
remains large at the solar age (see Fig. IT3| . This was shown independently by 
two groups, using slightly different approaches to rotational mixin g 11 ! and three 
different stellar evolution codes (Pinsonneault et al. 1989, Chaboyer et al. 1995, 
Talon 1997, Matias & Zahn 1998). 



10 The so-called lithium dip, discovered by Boesgaard & Tripicco (1986) and observed in all 
galactic clusters older than fs 300 Myr, refers to a strong lithium depletion in a narrow region 
of ~ 300 K in effective temperature, centered around T e g 6700 K. Lithium is burned at a 
temperature of ~ 2.5 X 10 6 K, which, in these stars, is located well into the radiative zone. As one 
looks at cooler stars, this burning layer is located closer to the base of the convective envelope. 

11 The Yale group uses a diffusive approximation for the meridional circulation instead of the 
formulation given in these notes. 
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Fig. 12. Predicted (closed circles) and observed (open circles; inverted open triangles 
are upper limits) lithium abundances in the Hyades around the lithium dip. Theoretical 
models are made with initial velocities of 50, 100 and 150 kms~ , and final velocities 
corresponding to those measured in the Hyades for the corresponding effective temper- 
ature. For cooler stars, lithium depletion remains large (JV(Li) is too small), at a level 
incompatible with observations. The inverted filled triangle corresponds to the prediction 
of rotational mixing if the inner rotation profile is artificially assumed to be solid body 
rotation, with surface braking. (Adapted from Fig. 6 of Talon & Charbonnel 1998.) 



These results on low-mass stars indicate that there must be another efficient 
transport process for angular momentum active when the surface convection zone 
appears; this would explain both the solar rotation profile and the rise of the 
lithium abundance on the cool side of the dip (Talon & Charbonnel 1998). 



2.6 Open problems in rotational mixing 

The 1-D modeling of meridional circulation is certainly a very important step in 
understanding and testing models of stellar evolution with rotation. However, one 
must not forget that it is actually an advective transport and that unidimensional 
modeling remains an approximation. 

The major limitation of the model comes from the invoked horizontal turbulent 
diffusion coefficient Dh- While in actual models it is linked to the source of hori- 
zontal shear (that is, some function of it) its magnitude remains quite uncertain. 
Furthermore, the hypothesis is that its effect will be to smooth out horizontal 
fluctuations. This is the case in the Earth's atmosphere albeit only up to scales 
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Fig. 13. Evolution of the solar rotation profile under the effect of rotational mixing only. 
The initial rotation velocity is assumed to be fOO kms -1 , and labels indicate the age in 
Gyr. (Figure 11.1 from Talon 1997.) 
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Fig. 14. Strongly anisotropic turbulence reduces the vertical extent of eddies, leading to 
a reduction of vertical transport. 



that are smaller than the Rossby radiua I. which is given by 



-^Rossby — 



H P N 
2D, sin 9 



(2.23) 



where 9 is the latitude. With typical values for Hp and N, it is of order Ljtossby — 
(50i?)/(u rot sin#), with v Iot in kms -1 . This limit will thus become relevant for 
stars rotating faster than about 100 kms -1 . 

Another issue concerns the vertical transport of chemicals in the presence of 
strongly anisotropic turbulence. As a turbulent eddy is displaced vertically in the 
fluid, it looses some of its coherence by being mixed by horizontal turbulence (see 
Fig. rrjj. This effect has been demonstrated in numerical simulations (Vincent 



This scale is the one at which Coriolis' force becomes comparable to vertical stratification. 
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et al. 1996, Toque et al. 2006). Talon, Richard & Michaud (2006) invoked such 
erosion in the context of AmFm stars. In § 11.51 we mentioned that AmFm stars can 
be understood as the result of radiative accelerations. However, there is a need for 
a competing process to reduce the magnitude of over/underabundances predicted 
in the case of a perfectly stable envelope. When rotational mixing prescriptions 
as described here are applied to such stars, the associated vertical transport of 
elements is too strong to account for AmFm stars; a reduction of the vertical 
transport by horizontal turbulence is required to correctly reproduce the observed 
abundances (Talon et al. 2006). However, such a reduction would degrade the good 
agreement obtained in the case of massive stars. Also, the effect of such anisotropic 
turbulence on a non-passive component (e. g. angular momentum) remains to be 
studied. 

A better understanding of the interaction between meridional circulation and 
turbulent properties could be obtained via the resolution of the full Navier-Stokes 
equations of hydrodynamics in 3-D (see also the 2-D, steady, calculations by Ga- 
raud 2002). This is a numerical challenge, as many time-scales are involved in 
the problem. A first step was taken in that direction by Talon et al. (2003) who 
wrote a code to tackle that problem. Used in direct numerical simulations with 
unrealistically large viscosities, it does lead to stationary solutions. A comparison 
of these solutions with those of analytical models should permit a verification of 
their validity in the limit of large viscosities. The goal is then to reduce the viscosi- 
ties and introduce a sub-scale turbulence model in order to represent unresolved 
features and to examine how solutions are modified. 



3 Magnetic fields 

3.1 Introduction 

There is a strong connection between rotation and magnetic fields. To examine 
that relation, let us start by taking a quick look at the equations that govern this 
interaction. We first have the induction equation 

— = V x {U x B- ryV x Bj , (3.1) 

which describes the action of a velocity field U on the magnetic field B in the 
presence of magnetic diffusivity -q. This equation is obtained by a combination of 
Ampere's law in the limit of no electric field and of Faraday's law of induction. We 
must also consider the momentum equation which takes into account the magnetic 
field by the means of the Lorentz force 
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The relative importance of magnetic field compared to inertial forces is given by 
the magnetic Reynolds number Re m 

ud 

Re m = — (3.3) 
V 

where u are d are respectively typical velocity and length scales. The "stretching" 
of field lines is efficient when this number is large. 

Here, we will discuss some of the impacts of magnetic fields on stellar structure, 
with special focus on their interaction with rotation. 



3.2 Ferraro's law of isorotation 



Probably the best known effect of magnetic fields is to produce solid body rotation: 
this is a restatement of Ferraro's theorem (1937), which we will examine here. To 
derive this theorem, we first separate our magnetic field into poloidal and toroidal 
components 



B = Br, 



Ba,€a, 



(3.4) 



with e0 a unit vector in the <j> direction. We now make various assumptions: 

• the fluid is perfectly conducting (r] = 0); 

• the fluid is incompressible (V ■ U = 0); 

• the solution is axisymmetric; 

• there is no meridional velocity field; 

• the solution is stationary. 

The induction equation for the Ba, component then reads 



dt 



V x 



(u x B) 

U (v • i?) - B (V • U) + (B ■ v) U - (u ■ v) B 
[B p ■ v) - (U^ ■ V) B 



(3.5) 



This equation shows how differential rotation produces a toroidal field from a 
purely poloidal field. In order to obtain a stationary state, we require 



(Bp ■ V) 



ft = 



(3.6) 



with Ua, = rsm8tt. This equation indicates that, under equilibrium conditions, 
fluid parcels that are located along the same poloidal field line must have the same 
angular velocity. 



3.3 Alfven speed 

It is not sufficient to know what the equilibrium state is, one must also evaluate 
the timescale required to reach this equilibrium. Here, we want to know how a 
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field line, that is stretched perpendicularly, responds. This gives rise to an Alfven 
wave, and the speed of this propagation is thus the Alfven speed v J^l 

v A = B (4Trp)- 1/2 (3.7) 

(Alfven 1942). In the star, shear will propagate over a stellar radius R in an Alfven 
time 

ta = — (3.8) 

VA 

which also defines the Alfven frequency 

= — = -=■• (3.9) 

TA H 

Let us mention that, in the absence of dissipation, this shear will simply bounce 
back and forth, and iso-rotation along poloidal field lines would not be obtained 
(see e.g. Mestel & Weiss 1987). Furthermore, in order to reach solid body rotation, 
diffusion is also required to "connect" the fluid velocity between various field lines. 
These two effects contribute to a significant increase of the timescale required to 
actually reach solid body rotation. 

3. 4 Solar spin-down 

We already mentioned in § 12.4.31 the problem of magnetic braking that is observed 
in all stars that have a "significant" surface convection zone, in practice all stars 
cooler than T e g ^ 7000 K. This can be understood in terms of mass loss in 
a magnetized wind, as originally suggested by Schatzman (1962) and Weber & 
Davis (1967): it is related to the fact that, when magnetic energy is larger than 
the wind's kinetic energy, matter is bound to co-rotate with the star even at large 
distances, leading to a large angular momentum loss even for a modest mass loss. 
An order of magnitude calculation is as follows. The surface solar dipolar field is 

B 0r w 1 G. (3.10) 

Magnetic flux conservation imposes 

r 2 B r = R%B 0r (3.11) 

while mass conservation leads to 

pu r r 2 = poUQ r RQ (3-12) 

where u r is the radial velocity in the wind for the mass loss. The comparison of 
the kinetic energy to the magnetic energy is given by the Alfvenic Mach number 

2 _ Anpul _ 1 Anplulr _ p A ,„..„, 

D r P ^Or P 



13 This velocity is obtained by coupling linearized versions of Eqs. 113.11 1 and 113.21 1. 
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which we use to define the Alfven density pa, corresponding to the density at which 
Ma = 1. If the mass loss velocity u 0r is known, one can evaluate this density. Using 
the values at the Earth orbit, ue — 400 kms" 1 and pe — 12 x 10 -24 gem" 3 leads 
to a density of 

Pa — 5 x 10~ 21 gem -3 . (3.14) 

Building a steady wind model similar to Parker's (1960), Weber & Davis (1967) 
concluded that the corresponding Alfven radius is 

r A =s 24i? . (3.15) 

The wind must co-rotate with the magnetic field up to that distance, explaining the 
efficiency of the surface braking. A more complete description of this process which 
also discusses the role of various field configurations is given by Kawaler (1988, 
see also references therein). This produces surface braking of the star. Through 
convective motions, it is expected that this braking affects the whole convection 
zone. The problem of the transport of angular momentum in the radiative zone 
remains. 

3.5 The solar rotation 

Another interesting application of Ferraro's law lies in the spin down of the solar 
radiative zone. As already mentioned (§ 12. 5p . helioseismology concludes that the 
the Sun's deep interior is rotating more or less rigidly, and rotational mixing alone 
is not able to provide an explanation of this feature. An obvious solution to this 
problem is to invoke a fossil magnetic field, as described at length by Mestel & 
Weiss (1987). Charbonneau & MacGregor (1993, see also Barnes, Charbonneau, 
& MacGregor 1999) studied this effect in a 2-D solar model, assuming a fixed 
poloidal configuration. They applied surface braking to a star originally rotating 
at 50 In these simulations, as the surface is spun-down, the poloidal field is 
sheared; this builds up a toroidal field that will then, through the Lorentz force, 
oppose the shear. This toroidal field grows until the torque exerted by magnetic 
stresses 

r = r 3 {BrB^) (3.16) 

is strong enough to oppose the wind-mediated torque. At this stage, the core and 
the envelope fully couple, and spin-down together, with the overall differential ro- 
tation remaining small. This occurs for all values of the poloidal field, although 
the time-scale for this re-coupling to occur depends on B r . Several configurations 
have been tested for the poloidal field, with the field lines anchored in the convec- 
tion zone (configurations Dl and D2 of their models), or not (configuration D4) 
or barely touching (configuration D3). 

Initial calculations were performed with the convection zone rotating as a solid 
body. When latitudinal differential rotation in the convection zone is present, 
due to Ferraro's law, it imprints itself in the interior (MacGregor & Charbon- 
neau 1999). Agreement with helioseismic inversion require that the field be totally 
contained in the radiative zone (configuration D4). 
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Fig. 15. Temporal evolution, in a 3-D MHD model of the solar radiative zone, of the 
angular velocity Q (color contours and white lines) and the mean axisymmetric poloidal 
field (superimposed black lines) , which initially is buried deeply in the interior (Fig. 4 of 
Brun & Zahn 2006). a-d) sequence spanning 4.7 Gyr (in solar equivalent units) of the 
dipolar magnetic field in the presence of rotation and shear. The field lines gradually 
connect with the imposed top shear, thereby enforcing differential rotation at all depths 
above latitudes greater than ~ 40°, at the (equivalent) age of the Sun. Note also in 
frames c-d the distortion near the poles of the lines of constant Q, which is associated 
with the Tayler instability (see § 13. 7[) . 

Recent 3-D MHD numerical simulations in which the fossil poloidal field is 
allowed to evolve have been performed to study this mechanism in a fully consistent 
way, using the ASH code (Brun & Zahn 2006 and references therein) . In all cases 
studied by these authors, the fossil field is found to diffuse outwards and inexorably 
connects with the outer convection zone, imprinting its differential rotation to 
the radiative interior. Contrary to the expectation of Gough & Mclntyre (1998), 
the meridional flow that develops between the base of the convection zone and 
the dipole poloidal field is not able to stop its progression. It remains to be 
seem whether things would change if lower (more realistic) values of the magnetic 
diffusivity, thermal diffusivity and viscosity were used (this is currently out of 
reach numerically). 

3.6 Magnetic stability 

Another issue that needs concern stellar astrophysicists regarding magnetic fields 
is that of the magnetic stability. This is well illustrated by the following quote 
from Braithwaite & Spruit (2004): 

"Although there have been educated guesses as to what shape a field 
in a stable equilibrium might have, all configurations studied with the 
analytic methods available so far were found to be unstable on the 
Alfven timescale; it has been impossible to prove stability for even a 
single field configuration." 
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However, while the presence of a large scale magnetic field in the solar radiative 
zone remains unproven, such fields are actually measured in several stars, as e.g. in 
Ap stars. Braithwaite & Spruit made numerical simulations of such stars, in order 
to understand what the structure of their field should look like. In their 3-D 
simulation, they model a non-rotating n — 3 polytrope, in which they introduced 
an initial random magnetic field. The total energy Eb in that initial field is 
assumed rather small, that is Eb ~ 0.01-Et, with Et the thermal energy of the 
star. 

They found that, after a few Alfven times, it was indeed possible to find a 
stable magnetic structure, which has the following properties: 

• the field is approximately axisymmetric; 

• the poloidal field has the same amplitude as the toroidal field; 

• magnetic helicity is conserved. 

The properties of their magnetic field are in close agreement with those actually 
measured in Ap stars, giving strength to the idea that it is of fossil rather than 
dynamo origin. 

3. 7 The Tayler-Spruit instability 

The issue of magnetic stability in rotating stars has been studied theoretically in 
great details by several authors, and an extensive review of these can be found in 
Tayler (1973) and Pitts & Tayler (1985). These linear stability study are performed 
using perturbations of the equations for the conservation of mass, momentum 
(Eq. 13. 2|) and energy and the magnetic induction equation (Eq. 13. ip and verifying 
whether a lower energy state exists. Actual calculations are time consuming, and 
have to be made for each magnetic configuration. The relative efficiency of these 
instabilities was reexamined by Spruit (1999), and he concluded that, the pinch- 
type instability, which is now more widely known as the Tayler-Spruit instability 
in the stellar community, was the first to set in. We will review this process here. 

Starting from a seed poloidal field B r , differential rotation winds it up and 
generates a mainly azimuthal field B^, 3> B r . The most unstable mode for this 
horizontal field is then a non-axisymmetric "kink" mode, in which the displace- 
ment is mainly horizontal (to minimize work against the density stratification) 
and incompressible (to minimize work against pressure) (see Fig. I16j) . Together, 
these two conditions imply 

lh/lr « ih/ir > 1 (3.17) 

where lh and l r are respectively the horizontal and vertical length scales and ^ 
and £ r are the horizontal and vertical displacements. Ignoring thermal diffusivity, 
the work against buoyancy is 

W = \p£N 2 (3.18) 
while the kinetic energy of the displacement is 



k = t P e^ A * i„eM 



(3.19) 
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Fig. 16. In the kink instability, toroidal field loops around the rotation axis move (mostly 
horizontally) relative to each other. Adapted from Tayler (1973). 



which yields the instability condition 

lr < ^ (3.20) 

assuming lh — r. This limits the vertical amplitude of the instability. However, 
magnetic diffusivity damps those small scales. The characteristic time-scale for 
magnetic diffusion td = l^/v has to be larger than the growth time of the instability 
which, in the limit of strong Coriolis force is of order 

a B =oj 2 A /fl, (3.21) 

for the instability to experience a net growth. One finds another condition for the 
vertical scale 

ll > (3.22) 
The combination of these two conditions yields a criterion for magnetic instability 

1/2 



This condition, valid in the limit of zero thermal diffusivity Kt, was originally 
derived by Acheson (1978). It would correspond to a region where N is dominated 
by mean molecular weight gradients. When Kt is taken into account (or in regions 
where TV is dominated by the thermal stratification), the required work against 
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gravity is reduced by Kx/rj, as in the case of hydrodynamic instabilities (§ 12.2.41 
and !2.3.2")) . The modified instability criterion then reads 

(Spruit 1999). For a sufficiently large magnetic field, there is a range of radial 
scales that are unstable with the largest unstable scale limited by stratification 
(Eq. I3.20|) and the smallest scale, by magnetic diffusivity (Eq. I3.22|) . These small 
scales produce turbulent magnetic diffusivity given by the equality in Eqs. (|3.23|l 
and (pT24]) 

— - m^) 4 (£) 2 - «- {vT&T <"»> 

in cases with (0) and without (1) thermal diffusivity. 

Spruit (2002) made further assumptions to produce a model that can be used 
in stellar evolution codes. Equating the growth-rate of (Eq. I3.2ip with the 
amplification time of the Tayler-Spruit instability 



N 

with q = 



<91nr 



(3.26) 



he obtains a relation between the Alfven frequency (or equivalently, the amplitude 
of the magnetic field) and the differential rotation 

~n =q N (3 - 27) 

valid in the absence of thermal diffusivity and 

when Kt cannot be neglected. 

Feed back to the rotation profile occurs through magnetic stresses. Still as- 
suming B^ ^> B r , one gets 

Ba, = (A-Kp) 112 ruj A and B r = ^Bj, m—B^ (3.29) 

lh r 

(l r is obtained by assuming the equality in condition 13. 20p . Magnetic stresses are 
then of the form 

1 /o \ 4 /o\ 1 ' 2 /fc \ 1 ' 2 

1 „2 n 2„3 il \ c 2 n 2 I S M I K T 



S = -B r B, = P r^ q ^-) or Sl = P r^ q ^-j (3.30) 
and this corresponds to a viscosity for the vertical transport of angular momentum 

s 2n2 /fi\ 4 __ 2r , fn\ 1/2 f k t \ 1/2 



^ = ^ = rfi<? [n) ot ™ = ' q {n K^J (3 ' 31) 
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Fig. 17. Rotation profile at the solar age for a model with meridional circulation and 
shear instability (dashed line) and with also the contribution of the Tayler-Spruit in- 
stability (dotted line). The initial velocity is 50 kms -1 , and surface braking follows 
Kawaler (1988). Dots with error bars correspond to the helioseismic solar rotation (Cou- 
vidat et al. 2003). From Eggenberger et al. (2005). 



with indexes (0) and (1) representing cases with and without thermal diffusivity 
(Spruit 2002). Maeder & Meynet (2004) made a more general formulation, which 
takes into account intermediate regimes. 

3.7.1 Applications to stars 

This model has been applied by the Geneva group to a variety of stellar models. 
This process is highly efficient in reducing differential rotation. In the case of mas- 
sive stars, rotation gradients are strongly reduced by the Tayler-Spruit instability 
(Maeder & Meynet 2004). Shear instabilities become much less important in this 
case, while meridional circulation is revived, since the rotation profile remains far 
from the one required by thermal equilibrium (§ !2.4.2p . These results are discussed 
at length by Meynet (2007, this volume). 

In the case of the Sun, the Tayler-Spruit instability also strongly reduces dif- 
ferential rotation. This was tested in an evolutionary model of a 1 Mq model, in 
which are also included meridional circulation and shear instability (Eggenberger, 
Maeder & Meynet 2005). This study confirms the capability of this process to 
reproduce the Sun's flat rotation profile (see Fig. |T7|) . 

3.8 Open problems in the coupling of magnetic fields and rotation 

The existence of the Tayler-Spruit instability has been observed in numerical sim- 
ulations (Braithwaite 2006). However, the interaction of this instability with the 
star's rotation with the mechanism described above implicitly assumes the regen- 
eration of a large scale, axisymmetric poloidal field from which the toroidal field 
can be generated. However, the "kink" instability described here corresponds to 
a global to = 1 mode, and so generates a poloidal field with a similar angular 
dependence, and not the axisymmetric field that is required in the model. This 
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was also shown in numerical simulations by Zahn, Mathis & Brun (2007). In their 
simulations, the dynamo for the poloidal field assumed by Spruit (2002) is not 
present. It remains to be seen if it would occur in conditions with smaller values 
of rj and v as expected in stellar interiors. 

Another approach to the interaction of magnetic fields and rotation consists 
in adding the Lorentz force in the equations for meridional circulation (Mathis & 
Zahn 2005). This leads to a complex set of 16 coupled differential equations that 
has not yet been included in actual stellar evolution calculations. 

4 Internal gravity waves 

4.1 Introduction 

A stratified fluid is one in which the density changes with height. In the oceans, 
such stratifications are stable if there is a continuous decrease in density with 
height; that is, lighter fluid is always located above heavier fluid. In the Earth's 
atmosphere and in stars, one has to remove the adiabatic temperature gradient 
when making this comparison. This gives rise to the well known Schwarzschild 
criterion. 

When a stable stratification is disturbed by a perturbation, buoyancy and 
gravity tend to restore the fluid to its previous equilibrium. For example, when 
a fluid parcel is displaced upward into a less dense region, gravity pulls it back 
down, below its equilibrium position and thus, into a denser region. Buoyancy 
takes over and pushes it upward. The fluid parcel "overshoots" to the less dense 
region and the cycle repeats. This oscillatory behavior forms what is called an 
internal gravity wave. 

In the Earth's atmosphere, wave induced momentum transport is a key element 
in the understanding of several phenomena; the most famous is the quasi-biennial 
oscillation (or QBO), which is an alternating pattern of westerly and easterly mean 
zonal winds observed in the stratosphere close to the equator. 

In astrophysics, internal gravity waves have initially been invoked as a source 
of mixing for elements induced by a variety of mechanisms (Press 1981, Garcia 
Lopez & Spruit 1991, Schatzman 1993, Montalban 1994, Montalban & Schatz- 
man 1996, 2000, Young et al. 2003, see § IO|l . Ando (1986) studied the transport 
of momentum associated with standing gravity waves; he showed how angular mo- 
mentum redistribution by these waves may increase the surface velocity to induce 
episodic mass- loss in Be stars (see also Lee 2006). He was the first to clearly state 
(in the stellar context) the fact that waves carry angular momentum from the 
region where they are excited to the region where they are dissipated. Goldreich 
& Nicholson (1989, see also Zahn 1975) later invoked gravity waves in order to 
explain the evolution of the velocity of binary stars, producing synchronization 
that proceeds from the surface to the core. Traveling internal gravity waves have 
since been invoked as an important process in the redistribution of angular mo- 
mentum in single stars spun down by a magnetic torque (Schatzman 1993, Kumar 
& Quataert 1997, Zahn, Talon, & Matias 1997, Talon, Kumar, & Zahn 2002). 
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Fig. 18. Numerical simulation of an internal gravity wave excited by an oscillatory forced 
patch. Colors show density fluctuations. Credit: Bruce R. Sutherland, U. of Alberta. 

4.2 Properties of internal gravity waves 

Internal gravity waves (IGWs) propagate due to buoyancy forces. Contrary to 
surface waves which are confined to a fluid's interface, internal gravity waves are 
able to propagate both horizontally and vertically through the fluid. They exist 
in stratified environments and transport energy and momentum. 

Because their restoring force is buoyancy, these waves are anisotropic with 
respect to "vertical" and "horizontal" orientations. This implies that, contrary to 
the case of a pressure wave, a point source does not produce an isotropic wave. 
Their frequency a has to be lower than the Brunt- Vaisala frequency N, which is 
the natural oscillation frequency of a displaced clement in a stratified region. We 
thus have 

0<a<N (4.1) 

(N is defined in § 12.2. ip . In stars, Nt is typically of order 10 _3 s -1 . A surprising 
property of internal waves is that their group velocity v g , which corresponds to 
the transport of energy is perpendicular to the phase velocity v p and to the wave 
numbei0 k. This is illustrated in Fig. [Tg] 

Let us now take a more formal look at the equations which describe those waves 
in a differentially rotating staiFI We use spherical coordinates (V, 9, cj>). e z is the 
unit vector along the rotation axis, and that in the azimuthal direction. We 



14 This can be demonstrated using the real part of the dispersion relation l|4.22|l . 
15 Our development follows Press (1981), Goldreich & Nicholson (1989) and Zahn et al. (1997). 
See also Unno et al. (1989) for details regarding the derivation of the wave equation. 
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assume that the angular velocity SI depends only on depth (just as in § I2.4|) . In 
an inertial frame, the total velocity U is 

U(r, 9, 0, t) = n(r) e z x r + u(r, 9, <j>, t), (4.2) 

where u is the velocity associated with the wave. The equation of motion reads 

^ + ^2Cle z x u + e^rsaxB (u ■ Vfi) j = --VP' + —g, (4.3) 



with 

d / d „ d 



at-\df n M)> (44) 

and where P' and p' are the Eulerian perturbations of pressure and density re- 
spectively. We used the Cowling (1941) approximation and neglected fluctuations 
of the gravitational potential. Viscosity is also ignored. We also consider the 
continuity equation 

^- + V • pu = 0, (4.5) 
at 

and the thermodynamic relation between p' and P' which in the adiabatic limit 
reads 

o' 1 P' N 2 

?--Tr-p-— C, = 0, (4.6) 
P Ti P g 

Ti being the adiabatic exponent and £ r the radial part of the displacement and 
which is related to the wave velocity by u r = ia^ r . 

We will proceed with a further simplification. Neglecting the terms in curly 
brackets in the equation of motion (|4.3[) , we treat waves as if they were pure 
gravity waves which are not modified by the Coriolis acceleration, but just feel the 
entrainment by the differential rotatiorF^I. In this case, we recover the equations 
that describe internal waves in a non-rotating star up to replacing the derivative 
with respect to time by (|4.4[) . see e.g. Unno et al. (1989). Solutions are then 
separable with respect to r, 9, (f> and t and are of the form 

u r (r, 9, 4>, t) = u v (r) P™ (cos 9) exp i (at - m<j>) , (4.7) 

where 

o-(r) = cr - mfl(r), (4.8) 

with <7o being the frequency in the inertial frame and P™ is the associated Leg- 
endre polynomial of order i and azimuthal number m. In the internal waves low 



1B The case of the gravito-inertial wave can be treated using the so-called traditional approx- 
imation which consists of neglecting the radial component of the Coriolis acceleration in the 
equation of motion. In that case, the horizontal structure can be described by Hough func- 
tions. This is discussed by Lee & Saio (1997) and Townsend (2003). The dispersion relation for 
gravito-inertial-Alfven waves is discussed in Kumar, Talon, & Zahn (1999). 
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frequency range which is considered here, terms of order a 2 can be neglected; this 
corresponds to the anelastic approximation. The function ^(r) = p^r 2 u v then 
obeys the second order equation 

d 2 * , ( N 2 ,\ £(£+1 



1 v ' y * = (4.9) 



dr 2 

and the dispersion relation is simply 

fc 2 cr 2 + N 2 kl = (4.10) 

(Unno et al. 1989). Here, we neglected a right hand side term of order H~ 2 which 
would mainly introduce phase shifts in the oscillation (cf. Press 1981). For future 
use, we introduce the vertical and horizontal wavenumbers, k v and kh, 



,., . N 2 \ t(l+l) ( N 2 . , 



This relation shows the highly anisotropic nature of internal waves, especially in 
the stellar context where we are mainly looking at waves with a <C N. When the 
wave is propagating, £ is conserved, kh varies as 1/r 2 and k v is modified according 
to (|4.11[) when the local frequency a and/or the Brunt- Vaisala frequency N change. 
If a > N, k v becomes imaginary and the wave is reflected. When a -C N, we have 
rk v ^> 1 and therefore the differential equation (|4.9j) may be solved by the WKB 
method, which yields 

u r = C r~'p _ 3 - l) P"\cos9) cos (at -m<f>- [ k v dr) , (4.12) 



with r c designating the base of the convection zonq^J and C is related to the 
wave amplitude. It describes a monochromatic wave of local frequency a, which 
propagate^ with a vertical group velocity 



dcr a N 2 - a 2 
dk~ v ~ ~~k v N 2 



(4.13) 



and horizontal velocities 



and 



rk v u r f . 

^ = - m ^TT)^ ( 4 - 14 ) 



'TV 2 

u 2 e +ul=[-Y-l) uf. (4.15) 



Here, we are implicitly assuming that IGWs are generated in the outer convection zone and 
travel inwards. 

18 The most general solution would include a stationary wave. 
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The <j) component of the velocity field depends on the azimuthal number m. In 
a rotating star, we distinguish between prograde (m > 0) and retrograde [m < 0) 
waves, corresponding respectively to waves traveling in the direction of rotation 
or against hi 19 l. 

So far, nonadiabatic effects have not been taken into account. To do so, we 
follow Press (1981) once more, and consider damping by thermal diffusion in the 
Boussinesq approximation^! and in Cartesian coordinate^]. In that case, the 
equation of motion reads 

i<rV 2 {pu r )= gkip' (4.16) 
while the equation of state may be written 



P_ 

P 



r 

) — 

T 



-ip- 



(4.17) 



where we again ignored the pressure perturbation. The fluctuation of fi is given 
by the conservation of mean molecular weight in the wave motion 



. A 4 V M 

1(7— = — — U r , 

fi Hp 



(4.18) 



and that of T by the heat equation 
V 



- Y = -(V ad -V)|- + fv 2 T', 



(4.19) 



where K is the thermal diffusivity. Combining Eqs. (|4.17|) . (|4.18p and (|4. 19[1 . we 
obtain 

K. 



N 2 

(ia - KV 2 ) p' = —pu r - — V" 
g to- 

and with (|4.16[) . we get the wave equation 



9 



pu r 



(4.20) 



,N 2 



K 



'4— j i>u, - ' — V J j V : ' -k 2 h ^r ) l'"r = 



g \ a 
From this equation, we derive the dispersion relation 



(4.21) 



k 2 - k 2 



1 



k 2 - p 



N 2 



(4.22) 



We now have Eq. (|4.9p which is the adiabatic equation in spherical geometry 
while Eq. (|4.2ip takes into account non-adiabaticity. We combine wanted features 



19 Note that some authors, including Unno et al. (1989), define the wave velocity as u r <x 
exp (im(f>) instead of u r <x exp (—im<j)) as used here. 

20 Here, the continuity equation l|4.5|l is replaced by V-u = 0, and implies that vertical variations 
of density are not considered in that equation. 

21 This derivation with mean molecular weight gradients is taken from Zahn et al. (1997). 



34 



Title : Stellar Nucleosynthesis: 50 years after B2FH 



from both these equations to obtain an equation that is locally Boussinesq and 
globally anelastic, and that takes into account nonadiabatic effects 

*$).-<>. <«»> 

As long as thermal damping is not too largj^l. the WKB solution for this equation 
is still given by equation (|4. 1 2[) . but multiplied by a damping factor exp(— r/2) 
where 



4.3 Excitation of internal waves 

In the Earth's atmosphere, IGWs are generated by several processes, like e.g. conv- 
ection and wind compression by topography. In the single star context, we are 
interested by waves produced by the injection of kinetic energy from a turbulent 
region to an adjacent stable regional. In the Earth's atmosphere, this is observed 
at the border of clouds (Townsend 1965). The existence of those waves in stars 
is observed in numerical simulations of penetrative convection both in 2-D and 
3-D (Hurlburt et al. 1986, 1994, Andersen 1994, Nordlund et al. 1996, Kiraga et 
al. 2000, Dintrans et al. 2005, Rogers & Glatzmeier 2005a, b). They are excited by 
two different processes namely 

• convective overshooting in a stable region; 

• excitation by the Reynolds stresses in the convection zone itself. 

Both sources contribute to the excitation. 

Ultimately, one would like to obtain realistic wave fluxes from numerical sim- 
ulation themselves. Some improvement for such prescriptions has been made by 
Rogers & Glatzmeier (2005b) who calculated wave excitation in a cylindrical (2- 
D) model having a stratification that is similar to that of the Sun. However, the 
level of turbulence reached in that simulation is far from realistic and thus, cannot 
provide a good answer to the question yet. Furthermore, an analysis of the wave 
spectrum characteristics as a function of convective properties (such as the con- 
vective luminosity and/or the turn-over time-scale) has not been done. For stellar 
applications, the variation of wave generation with effective temperature^ is a key 
feature in understanding their differential properties (see § I4.6.2| . These are the 
two reasons why we will rather rely on theoretical estimates for wave generation. 



If the imaginary part of this equation is larger than the real part, the solution for u r is 
evanescent rather than oscillatory. 

23 In binary stars, the variation of the gravitational potential can excite a gravity wave in early 
type stars. Its effect has been studied by Zahn (1975a) and Goldreich & Nicholson (1989). 

24 One should note that the properties of the surface convection zone, like its depth and relevant 
time-scales, are determined almost uniquely by the star's surface temperature. 
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Excitation by overshooting is quite difficult to evaluate analytically. Garcia 
Lopez & Spruit (1991) made such an attempt and assumed that the pressure 
perturbation produced by turbulent eddies at the radiative/convective boundary 
is equal to the wave pressure perturbation. In this model, small stochastic eddies 
combine to excite a whole wavelength spectrum. It has been formulated under the 
assumption of homogeneous turbulence, but it could be modified to incorporate 
the presence of down-drafts as observed in numerical simulations. Comparisons 
of numerical results with this theoretical model show that the maximum in the 
wave spectrum produced by the simulation was similar in amplitude to that of 
this parametric model (Kiraga et al. 2003, Rogers & Glatzmeier 2005b). However, 
in simulations, modes were excited over a much broader range of frequencies and 
wavelengths. It is not clear whether this is related to the bi-dimensional nature of 
the simulations or to shortcomings in the model. 

Another estimate for excitation by overshooting has been made by Fritts, 
Vadas, & Andreassen (1998), in a study aimed at estimating the residual circula- 
tion induced by latitude dependent wave dissipation in the tachocline. However, 
they concentrate on the small wavelength waves that dissipate close to the con- 
vection zone, and their mechanism does not take the combination of small scale 
eddies into account in order to produce low degree waves. Such waves are essential 
in order to influence the inner regions on an evolutionary time-scale (§ I4.5( see also 
Talon, Kumar, & Zahn 2002). Let us also note that such low degree waves are 
actually observed in the wave spectra of numerical simulations. 

IGWs can also be excited in the convection zone itself. That can be estimated 
by convolving the local wave function, which is evanescent and whose local am- 
plitude is proportional to exp [— J dr fc r ] , with entropy fluctuations and Reynolds 
stresses. A complete description of this process has been developed by Goldreich 
and his collaborators (Goldreich & Keeley 1977, Goldreich & Kumar 1990, Gol- 
dreich, Murray, & Kumar 1994, hereafter GMK). Their model was first applied 
to the solar p-modes, and quite successfully reproduces the solar spectral energy 
input rate distribution, provided one free parameter which describes the geometry 
of turbulent eddies is calibrated. In that case, driving is dominated by entropy 
fluctuations. Balmforth (1992) made a similar study, using a somewhat different 
formalism. Subject to the calibration of a free parameter, he is also able to re- 
produce the spectral energy distribution; however, Reynolds-stresses dominate the 
driving in his calculations. 

The GMK formalism can be adapted to the case of g-modes (Kumar & Quataert 
1997); the energy flux per unit frequency Te is then 



x exp 



i) 



2r 2 



3r4 



l + (ar L y 



(4.25) 



where £ r and [^(£+1)] 2 £ h are the radial and horizontal displacement wave- functions, 
which are normalized to unit energy flux just below the convection zone, v is the 
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Fig. 19. Penetration of IGW into the convection zone of a ZAMS 1 Mq model for various 
degrees I. The arrows indicate the depth where the amplitude is reduced by a factor of 
2. Also shown is the local convective flux given by the mixing length theory. 



convective velocity, L the radial size of an energy bearing turbulent eddy, T£ « L/v 
the characteristic convective time, and h a is the radial size of the largest eddy at 
depth r with characteristic frequency of a or greater (h a — Lmin{l, (2<ttx)~ J } ). 

In the convection zone, the mode is evanescent and the penetration deptH 25 ! is 
inversely proportional to y/£(£ + 1). Figure [19] compares this distance for various 
modes with the convective energy that is locally available. In the outer part of the 
model, the luminosity is carried almost uniquely by convection. 

4.4 Particle transport 

While we all learned in basic physics lectures that waves do not produce transport 
of matter, this is actually only true when a first order, linear approximation is 
used. In fact, there exists a wide range of physical processes that may give rise 
to mixing by waves, and we will briefly review some of these here. Let us note 
however that the actual calculations that were made with these in a stellar context 
do not take into account angular momentum transport by the same waves, which 
would greatly reduce their penetration depth (see § 14. 5p and thus, the (numerical) 
results that are mentioned should be reviewed accordingly. 

4.4.1 Shear unstable waves 

Let us look at the case of low frequency IGWs in which a <C N. Then, the 
velocity of the associated flow is nearly horizontal, and produces a sinusoidally 
varying vertical shear. Such a flow has inflection points, which may trigger a shear 



25 This is the theoretical dependence. Damping can be enhanced by turbulence, which would 
reduce the surface amplitude of the mode (see e.g. Andersen 1996). 
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ABC 




Fig. 20. Movement of a fluid parcel (whose initial position is shown with a gray dot) 
when a wave passes in the linear regime (A), when thermal diffusion is taken into account 
(B) and when non-linearities are taken into account (C). In (B), the net displacement is 
in the direction of gravity while in (C), it is perpendicular to g. 



instability that will produce small scale turbulence. Because the flow is time- 
dependent, the shear rate k v Uh has to be larger than the wave frequency for this 
to occur and hence the condition for weak mixing by this process is (Press 1981, 
Garcia Lopez & Spruit 1991) 

k v u h > o. (4.26) 

When this condition is satisfied, the local shear rate can be compared to the sta- 
bilizing Brunt- Valsala frequency, leading to the well known Richardson criterion. 
Instability occurs when 

N 2 1 

{k v u h f < 4 



Ri = ^ < -■ ( 4 - 27 ) 



Taking thermal diffusivity into account, the Richardson criterion is modified and 
becomes 

(k v Uh) K 4 

where iovo = Do is the turbulent diffusion associated with eddies of size £o and 
velocity vo- The actual turbulent diffusion D wavc s hoar will be given by the largest 
value of this product that is unstable and thus 

^wavc shear — ^ ^2 (4.29) 

Garcia Lopez & Spruit (1991) applied such a formula to the case of F stars, 
in order to reproduce the so-called Li dip (see § I4.6.2|) . Using an excitation model 
based on the matching of pressure fluctuations at the base of the convection zone 
and that of the waves, they showed that Li burning could be produced in the dip 
region if the input of turbulent energy in waves is large enough 26 ! and provided the 
mixing length of convection is suitably chosen. 



26 In their paper, Garcia Lopez & Spruit multiplied the wave flux calculated using a standard 
mixing length approach of convection by 10 to obtain enough depletion in the Li-dip. This could 
be explained by considering the presence of plumes in convection. 
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4.4.2 Irreversible second-order stochastic diffusion 

In the adiabatic and linear limit, the trajectory of a fluid parcel as a wave passes is 
a simple circular movement (Fig. I20I A.V However, when thermal diffusion is taken 
into account, the "memory" of the initial entropy content is lost, and the fluid 
parcel does not return to its initial entropy level (Fig. I20B: Press 1981, Schatz- 
man 1993, Montalban 1994). In the adiabatic limit, the vertical displacement of 
the fluid element is given by 6 V — u v /a, while in the absence of a restoring force 
it would be A v = k v 8 2 . The fraction of entropy that is lost during a vertical 
excursion can be roughly estimated as the fraction of this length-scale A„ corre- 
sponding to the ratio of the diffusive time-scale K k 2 , and the period 1/rr. Overall, 
the vertical displacement of the fluid parcel is 

ju3 ^ .2 

I = (4.30) 

(Press 1981). The velocity scale associated with this process is 

v = ta, (4.31) 
and this leads to a diffusion coefficient of the form 

D 2nd order = Iv = K 2 ^ (4.32) 

for one wave. A combination of all excited waves has to be made to transform this 
into a global diffusion coefficient (this is discussed at length by Montalban 1994). 
This process should be efficient only over about one damping length scale and thus, 
will lead to a diffusion coefficient that decreases quite rapidly with depth. Calcu- 
lations using such a description were made by Montalban & Schatzman (1996) to 
compute Li destruction in cool main sequence stars. 



4.4.3 Stokes' drift 

Stokes ' drift can be defined as the difference between the mean fluid velocity mea- 
sured at a point and the mean velocity of a drifter in the direction of wave prop- 
agation and was originally described by Stokes (1847, see Lighthill 1979 for a 
complete description). If there is no mean velocity, it represents the velocity of 
a passive scalar that is slowly moving through the fluid. This displacement is in 
the direction of wave propagation (v^>), which is mainly horizontal for the IGWs 
we are studying here. Heuristically, this drift is caused by an asymmetry between 
the horizontal displacement when the wave is moving downward and upward. The 
resulting movement is shown in Fig. [20lC. In the case of stars, this will produce 
mainly a horizontal redistribution of elements and thus, little vertical mixing is 
expected from this process. 
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4.4.4 Wave-induced turbulence 

The last process I will describe here is wave-induced turbulence. This has been 
discussed by a few authors, under various forms. Canuto (2002) developed a model 
in which gravity waves act as a source term in the equation that describes turbu- 
lence. Young et al. (2003) describe mixing by waves generated at the boundary of 
a convective region as a physical description of "overshooting" . Talon & Charbon- 
nel (2005) study the formation of an oscillating shear layer caused by differential 
wave damping (see § l4.5.ip . In that process, energy is transfered by waves from 
the convection zone to the shear layer and stored in differential rotation, which 
can then be converted to turbulence by the shear instability (§ I2.3.2|) . Let us note 
that it is different from the process described in § 14.4.11 here it is the rotation 
profile that becomes unstable, and not the wave itself. 

In all these cases, turbulent mixing should be large in a rather small region 
just below the convection zone, and acts to smooth out the possible chemical 
discontinuity at convective boundaries. 



4.5 Momentum transport 

Let us now describe how angular momentum is transported by IGwJ^I. First, the 
horizontal average of the kinetic energy density can be calculated from the WKB 
velocities (|4.12l) and (|4.15j) and yield 

^p<u 2 > = -L JJ ^ p( u 2 r + u 2 g + ul) sin9 d9d(f> 

1 N 2 1 f f IN 2 
= ~- I —JJpu 2 r ^ddedcjy=-—p<u 2 r >. (4.33) 

We now multiply this quantity by the group velocity (|4.13j) to get the average 
kinetic energy flux carried by a traveling IGW 

T K ^, m = - P < u > Yv ^^ = g P < u > W2 . (4.34) 

The mean radial flux of angular momentum is evaluated using Eqs. (|4.12p . 
flOp) , fi^33l) and (14341) 

3~j.a,i,m = — J J P t sin 6 UtjjUr sin d0 d(j) 

= ~ m r JTTV TT J" // P U r Sm 9 d9 ^ = 2 - ^K.a.Lm- (4.35) 
+ 4tt J J (7 

Inserting (|4. 1 2[) into (|4.35[) . it is possible to verify that, in the adiabatic limit, the 



27 This description follows Zahn et al. (1997). 
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Fig. 21. Angular momentum wave luminosity C.j, a ,i in a ZAMS IMq model for various 
degrees £ and frequencies v = <j/2tt according to the GMK model. 



angular momentum luminosity is conserve* 



,28 



est. 



(4.36) 



Figure [21] shows the wave spectrum of angular momentum luminosity in a 1 Mq 
model when the GMK model is used for the kinetic energy flux. 

When radiative damping is taken into account, the wave amplitude is multi- 
plied by an attenuation factor exp(— r/2) (cf. Eq. [4.24(1 . Locally, the total angular 
momentum luminosity is the sum of the contribution of all waves 



£.i( r ) = ^2 £j,v,£,m( r c) exp [-r(r, a, £)] 



(4.37) 



The deposition of angular momentum is then given by the radial derivative of this 
luminosity. Let us remark that here we consider only the radial dependency of 
wave transport, thus the use of horizontal averaging. When other processes for 
angular momentum redistribution are also taken into account, the evolution of 
angular momentum follows 



dt 



J_d_ 

5r 2 dr 



[ P r 4 nu] 



J_d_ 

r 2 dr 



pv t r 



dr 



± ,4 - 38) 



where U is the radial meridional circulation velocity and v t the viscosity from 
hydrodynamical instabilities (§ I2.2|f^l . The "+" ("— ") sign in front of the angular 
momentum luminosity corresponds to a wave traveling inward (outward). 



8 On the other hand, the kinetic energy of the wave varies with depth since the local frequency 
a depends on the rotation rate Q(r). Bretherton (1969) showed that, in the plane-parallel case, 
conservation is ensured in the inertial frame. 

29 One could also add the magnetic viscosity of § 13.71 but for consistency gravito- Alfven waves 
should then be considered. See e.g. Barnes, MacGregor & Charbonneau 1998. 
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Fig. 22. Evaluation of the damping factor r (cf. Eq |4.24[) for a frequency of 1 fiHz in 
a solar model. The depth corresponding to an attenuation by a factor 1/e is shown for 
various degrees I. 
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Fig. 23. Shear layer oscillation (SLO) in a 1.3 Mq model. Successive profiles are sep- 
arated by 1 year intervals. The layer's thickness depends on the magnitude of thermal 
dissipation and the oscillation period is a function of the total wave angular momentum 
luminosity. 



4.5.1 Wave-mean flow interaction and the SLO 

We now examine the impact of IGWs on the radial distribution of angular mo- 
mentum. Figure [22] shows the penetration depth for waves of various degrees. For 
low frequency waves, damping is strong enough so that most waves are damped 
very close to the base of the convection zone. These are the waves we will examine 
here. 

Let us first take a look at the damping integral (|4.24[) . and assume that both 
prograde and retrograde waves are excited with the same amplitude (according to 
Eq. l4.25|) . In solid body rotation, both waves are equally dissipated when traveling 
inward and there is no impact on the distribution of angular momentum. In the 
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presence of differential rotation, the situation is different. If the interior is rotating 
faster than the convection zone, the local frequency of prograde waves diminishes, 
which enhances their dissipation; the corresponding retrograde waves are then dis- 
sipated further inside. This produces an increase of the local differential rotation, 
and creates a double-peaked shear layei0. In the presence of shear turbulence, 
this layer oscillates, producing a "shear layer oscillation" or SLO (cf. Fig. [23l and 
Gough & Mclntyre 1998, Ringot 1998, Kumar, Talon, & Zahn 1999). Viscosity is 
essential in this process and is responsible for the disappearance of the outermost 
peak as its slope gets steeper and the peak is slowly absorbed by the convection 
zone. This oscillation is similar to the Earth's QBO (see e.g. Mclntyre 2003 for 
details regarding the QBO). 

4.5.2 Secular effects 

The SLO, by locally Doppler shifting to lower frequencies both prograde and retro- 
grade waves, acts as an efficient IGW filter (see Eq. I4.24|) and for some time there 
was doubt that any wave power would remain beyond (Gough & Mclntyre 1998). 
However, Talon, Kumar, & Zahn (2002) and Talon & Charbonnel (2005) per- 
formed complete calculations of wave dissipation in the presence of a SLO and 
showed that the low-degree waves do conserve a significant amplitude beyond the 
shear layer. 

If the radiative zone has the same rotation rate as the convection zone, over 
a complete SLO cycle, the magnitude of the prograde and retrograde peakqfj are 
on average equal, and waves of azimuthal number +m and —m have (on average) 
equal amplitudes after crossing the SLO. In the presence of differential rotation 
however, with the inner part rotating faster than the outer part as is the case 
of a solar type star that is spun down by magnetic torque, the prograde peak is 
always larger than the retrograde peak and there is a net flux of negative angular 
momentum to be redistributed in the star's radiative zone. Talon, Kumar, & 
Zahn (2002) showed that this can spin down the interior of a star over long time- 
scales. Talon & Charbonnel (2005) later showed that this asymmetric filtering 
depends only on the difference of rotation rates at the base of the convection zone 
and at the base of the SLO (5il — £l cz — S^slo)- 

Figure [24] shows how wave amplitudes are differentially reduced by crossing 
the shear layer. The rightmost column shows the net luminosity beyond the SLO; 
the low-frequency, low-degree waves have the strongest impact on the secular re- 
distribution of angular momentum in the radiative part of the star. They must be 
low-frequency, otherwise differential damping is too weak, and low-degree, other- 
wise they are damped in the SLO rather than beyond. 



30 Because local shears are amplified by waves, even a small perturbation will trigger this. 

31 The prograde peak refers to the layer that is rotating more rapidly than the convection zone 
and it produces filtering of prograde waves, while the retrograde peak is rotating slower and 
filters retrograde waves. 
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Fig. 24. Wave characteristics below the shear layer for a differential rotation of SQ = 
0.1 /iHz over 5% in radius in a 1.3 M star. First and second columns: Wave amplitude 
for m = ±1 (top) and m = ±£ (bottom). Third column: Amplitude difference between 
retrograde and prograde waves of the same azimuthal number. Fourth column: Net 
luminosity for a given azimuthal number. 
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Fig. 25. Evolution of the interior rotation profile in a 1 Mq model with meridional circu- 
lation, shear instabilities and IGWs. The initial equatorial rotation velocity is 50 kms -1 
and surface magnetic braking is applied. Curves correspond to ages of 0.2, 0.21, 0.22, 
0.23, 0.25,0.27 (dashed lines), 0.5, 0.7, 1.0, 1.5, 3.0 and 4.6 (solid lines) Gyr. Adapted 
from Charbonnel & Talon (2005). 
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Fig. 26. Evolution of the surface lithium abundance for various initial rotation velocities 
in models including meridional circulation, shear instabilities and IGWs. The surface Li 
abundances as measured in various clusters are also shown. The vertical extent of boxes 
shows the range of lithium values for stars with an effective temperature correspond- 
ing to that of the model ±100 K at the cluster age plus a typical error in abundance 
determination. Adapted from Charbonnel & Talon (2005). 

4.6 Applications of angular momentum transport 
4.6.1 The solar rotation 

This description of IGWs has been applied to the solar rotation problem. Talon, 
Kumar, & Zahn (2002) proved, in a static solar model, that IGWs could be respon- 
sible for the required angular momentum extraction. More recently, Charbonnel & 
Talon (2005) revisited the problem and looked at the role of IGWs in an evolving 
solar model, which undergoes spin-down from the zero age main sequence (ZAMS) 
on. They showed that, at the age of the Sun, differential rotation indeed becomes 
quite small (see Fig. |25|) . In this simulation, the low-degree waves penetrate all 
the way to the core and spin it down extremely efficiently at the very beginning of 
the evolution. This is due to the small amount of angular momentum (cx r 2 ) con- 
tained in the core. Once the core has been spun down, the damping of retrograde 
waves, that carry the negative angular momentum, increases locally. Consequently 
a "slowness" front forms and propagates in a wave-like way from the core to the 
surface. As further braking takes place, a second front forms and propagates out- 
ward. The differential wave filtering adjusts itself so as to compensate for the 
flux of angular momentum that is lost through the stellar wind. This explains 
why front propagation is fast at the beginning and then slows down, just as the 
spin-down rate does. 

Figure UHl shows the evolution of the surface Li abundance in a solar model 
for various initial velocities. One should note here the (small) dispersion that is 
obtained when very different initial velocities are considered. This would not be 
the case if only meridional circulation and shear turbulence where considered since 
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Fig. 27. (top) Lithium dip as observed in the Hyades. (bottom) Net angular momentum 
wave luminosity below the SLO for a differential rotation SQ = 0.05 /xHz over 0.05 R*. 
The vertical lines correspond to the requirements in the efficiency of angular momentum 
transport by a "non-standard" process so that rotational mixing would cause the Li dip 
as suggested by Talon & Charbonnel (1998). Adapted from Talon & Charbonnel 2003. 



then, mixing would be proportional to the angular momentum loss. 
4.6.2 The Li dip 

We mentioned (§ I2.5|) that rotational mixing could explain the hot side of the Li 
dip but destroy too much lithium when looking at cooler stars. TC98 suggested 
that if the process that is required to explain the Sun's rotation profile starts to 
become efficient in the center of the Li dip, these observations would be reconciled 
with rotational mixing, thus giving a model that is valid for all main sequence 
stars. One then identifies three regimes: 

• Stars with T e g Si 6900 K have a very shallow convective envelope which is 
not an efficient site for magnetic generation via a dynamo process. They are 
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thus not slowed down by a magnetic torque. As a result, these stars soon 
reach the equilibrium meridional circulation (§ |2.4.2[1 and the associated weak 
mixing is just sufficient to counteract atomic diffusion. 

• In stars with ~ 6900 ^ T e g ^ 6600 K, the convective envelope deepens and 
a weak magnetic torque appears that slows down the outer layers. In this 
case the transport of angular momentum by meridional circulation and shear 
turbulence increases, leading to a larger destruction of Li. 

• Stars on the cool side of the Li dip (T e g ^ 6600 K) have an even deeper 
convective envelope sustaining a very efficient dynamo which produces a 
strong magnetic torque that spins down the outer layers very efficiently. This 
outer convection zone is an efficient source for IGW generation and wave- 
induced transport reduces the efficiency of rotational mixing, thus producing 
a rise in Li abundances on the cold side of the Li dip. 

Talon & Charbonnel (2003) showed that the net angular momentum luminosity 
corresponding to IGWs does have the proper temperature dependency (Fig. |27|) 
and thus could well explain the cool side of the lithium dip. Complete evolutionary 
calculations are underway to verify this proposition. 

4. 7 Open problems in internal waves physics 

While results obtained including internal waves in evolutionary models are quite 
promising, they remain crude, and several points are yet to be addressed. First 
of all, one should revisit the excitation of IGWs. Recent progress has been made 
in the description of the excitation of solar p-modes by considering better models 
for convection (see e.g. Belkacem et al. 2006). Work is underway to adapt this 
description to IGWs (Belkacem et al. in prep.) Another important point is the 
inclusion of the Coriolis force in calculations. This issue has been raised by several 
authors (Talon 1997, Lee & Saio 1997, Townsend 2003 and Mathis 2005) and work 
is underway to correctly incorporate it in descriptions of momentum transport in 
massive stars (Pantillon, Talon & Charbonnel in prep.) Finally, all calculations 
that are presented here correspond to horizontal averages. Implicitly, this assumes 
that horizontal transport is efficient in the redistribution of angular momentum 
along isobars. This faces the "Rossby scale problem" already mentioned in § 12.61 
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